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Objectives: 

The  properties  of  high  energy  density  matter  (HEDM)  were  under  investigation  in  our  group 
Findings  are  described  below. 

A.  Low  Temperature  HEDM  Simulations 

Solid  hy(togen  doped  with  energetic  impurities  may  form  the  basis  for  HEDM  to  be  used  in 
rocket  propulsion.  A  key  research  priority  has  been  to  investigate  the  equilibrium  properties  of  a 
boron  impurity  and  of  an  aluminum  impurity  trapped  in  solid  /?ara-hydrogen.  The  simulation 
methods  used  in  the  study  have  allowed  us  to  take  into  account  the  quantum  nature  of  the  HEDM  in 
question.  We  strongly  believed  that  using  classical  molecular  dynamics  simulations  to  evaluate 
these  equilibrium  properties  would  predict  qualitatively  incorrect  results  for  low  temperature  solid 
hydrogen,  because  of  the  highly  quantum  nature  of  the  hydrogen  matrix. 

The  goal  of  our  group's  Air  Force-supported  low  temperature  HEDM  research  program  was 
to  understand  from  large-scale  computer  simulation  the  differences  in  the  equilibrium  properties  of 
the  HEDM  when  the  doped  impurity  is  aluminum,  and  to  asses  the  effects  of  anisotropic 
interactions  between  the  impurity  and  the  molecule.  In  low  temperature  hydrogen  HEDM, 
quantum  effects  are  enormous,  so  a  highly  specialized  method  was  required  for  these  demanding 
computer  simulations. 

The  primary  effort  over  the  course  of  this  project  was  the  large-scale  simulation  of  aluminum 
and  boron  impurities  in  para-hydrogtn  at  4K.  Significant  results  were  obtained  which  are 
summarized  in  the  next  section. 


Accomplishments/New  Findings: 

A.  Low  Temperature  HEDM  Simulations 

As  stated  earlier,  the  computer  simulation  of  impurity  trapped  in  solid  hydrogen  must  account 
for  the  highly  quantum  nature  of  the  system.  The  method  used  in  our  group  to  study  this  problem 
is  path  integral  molecular  dynamics  (PIMD).  In  PIMD  the  ensemble  of  a  given  number  of 
quantum  particles  is  represented  quantum  mechanically  by  a  certain  set  of  that  number  of 
pseudoparticles.  The  dynamics  were  performed  on  these  pseudoparticles  and  the  phase  space  was 
properly  sampled.  This  allowed  us  to  obtain  correct  results  for  the  equilibrium  properties  of  a 
highly  quantum  system  like  the  HEDM  studied  here. 

As  stated  above,  an  important  part  of  this  study  was  assessing  the  effect  of  the  anisotropic 
interactions  of  the  impurity  with  the  pHj  molecules  in  the  matrix.  This  directionality  in  the 
interaction  is  due  to  the  fact  that  aluminum  has  a  singly  filled  3p  orbital  in  its  ground  state. 

The  potential  of  the  Al-pHj  pair  is  then  described  by  two  degenerate  H  states  and  a  S"''  state. 
If  the  z-axis  is  taken  to  be  the  intemuclear  axis,  then  these  three  states  correspond  to  the  three 
Cartesian  orientations  of  the  singly  filled  p-orbital  {p^  for  the  2^  state  and  p^  and  p^  for  the  H  states). 

We  have  compared  the  calculated  equilibrium  properties  of  the  system  with  explicit  inclusion 
of  the  electronic  anisotropy  of  the  impurity-pHj  interaction,  and  with  a  simplified  potential  in  which 
the  impurity-pH2  pair  potential  is  averaged  over  the  three  possible  orientations  of  the  /^-orbital 
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Insight  into  the  structure  of  the  equilibrated  systems  is  given  by  the  impurity-pHj  pair 
distribution  functions,  g(R).  These  were  calculated  for  the  single  site  substituted  system  and  for  the 
system  containing  a  defect  in  the  nearest  neighbour  site  of  the  impurity.  The  site  substituted 
system  showed  less  structural  change  than  the  system  with  a  defect.  Shown  in  Figure  1  are  the 
g(R)s  for  the  case  of  a  aluminium  impurity  without  a  defect. 


Figure  1 

The  shift  in  the  profile  of  the  g(R)  and  the  change  in  the  structure  of  the  first  peak  in  Figure  1 
implies  that  there  is  some  structural  rearrangement  of  the  lattice  in  the  vicinity  of  the  impurity, 
and/or  rearrangement  of  the  impurity,  when  the  orientation  dependent  potential  is  used. 


iniportant  part  of  this  project  was  assessing  the  effect  of  the  anisotropic  interactions  of  the 
impurity  with  the  pHj  molecules  in  the  matrix.  This  directionality  in  the  interaction  is  due  to  the 
fact  that  both  boron  and  aluminum  have  a  singly  filled  p  orbital  (2p  in  the  case  of  B  and  in  the 
case  of  Al). 

For  exain^ple,  in  the  case  of  the  boron,  the  pairwise  interaction  of  the  B  in  its  2s^2p  electronic 
ground  state  with  the  spherical  partner  can  be  described  by  two  degenerate  IT  states  and  a  2;"  state. 
If  the  z-axis  is  t^en  to  be  the  intemuclear  axis,  then  these  three  states  correspond  to  the  three 
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For  each  case,  we  have  compared  the  calculated  equilibrium  properties  of  the  system  with 
explicit  inclusion  of  the  electronic  anisotropy  of  the  impurity-pHj  interaction,  and  with  a  simplified 
potential  in  which  the  impimty-pHj  pair  potential  is  averaged  over  the  three  possible  orientations  of 
the/?-orbital.  Shown  in  Figure  2  are  the  g(R)s  for  the  case  of  a  boron  impurity  with  a  defect. 


2 


3 


wo 


2 


1 


0 


r(A) 


Figure  2 

The  shift  in  the  profile  of  the  g(R)  in  Figure  2  implies  that  there  is  some  structural  rearrangement  of 
the  lattice  in  the  vicinity  of  the  impurity,  and/or  rearrangement  of  the  impurity,  when  the  orientation 
dependent  potential  is  used. 


The  effect  was  much  more  pronounced  in  the  case  of  an  aluminium  impurity  with  a  neighboring 
lattice  vacancy.  From  the  Al-pHj  g(R)  in  Figure  3,  it  was  apparent  that  a  significant  structural 
rearrangement  occured  when  the  orientation  dependent  potential  was  used.  The  two  peaks  seen  in 
the  profile  for  the  spherically  averaged  potential  split  into  two  when  the  orientation  dependent 
potential  was  used.  TOs  indicates  that  the  arrangement  of  the  molecules  around  the  A1  is  no 
longer  as  uniform  as  in  the  case  of  the  spherical  potential.  To  gain  insight  into  the  cause  of  the 
ch^ge  in  the  g(R)  profiles  discussed  above,  we  projected  onto  the  hexagonal  plane  of  the  hep 
lattice  the  positions  of  the  A1  atom  and  the  1 2  nearest  neighbor  pHj  molecules  for  the  case  of  the 
aluminum  impurity  with  an  adjacent  defect. 

In  Ae  site-substituted  hep  lattice,  six  of  the  nearest  neighbors  are  situated  at  the  vertices  of  a 
hexagon  in  the  xy  plane,  three  additional  pHj  molecules  are  located  at  the  vertices  of  an  equilateral 
triangle  above  the  xy  plane,  and  the  final  three  pHj  molecules  are  at  the  vertices  of  an  equilateral 
triangle  below  the  xy  plane.  Presented  in  figure  4  is  projection  obtained  using  the  spherically 
averaged  potential.  The  origin  is  defined  by  the  initial  hep  position  of  the  A1  atom.  The  filled 
circles  and  the  filled  tri^gle  are  the  positions  of  the  ideal  hep  lattice  sites.  The  squares  are  the 
averaged  centroid  positions  of  the  pHj  molecules.  The  open  triangle  is  the  averaged  centroid 
position  of  the  A1  atom. 
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Figure  4 

In  this  case  some  of  the  nearest  neighbor  hydrogen  molecules  move  from  their  hep  positions,  where 
some  stay  reasonably  close  to  them.  The  mean  A1  atom  position  shifts,  however,  slightly  towards 
the  vacancy.  The  projection  presented  in  figure  5  is  for  Ae  case  where  the  orientation  dependent 
potential  was  used.  In  this  plot  the  open  circle  represents  the  location  of  the  vacancy.  This 
projection  shows  that  upon  equilibration  the  A1  ends  up  in  between  its  substitution  site  and  the 
vacancy.  There  is  also  evidence  that  the  molecules  shift  from  their  hep  positions  slightly 


4 


towards  the  new  position  of  the  A1  atom.  This  large  shift  in  position  of  the  A1  impurity  arises  from 
the  strongly  anisotropic  interaction  with  the  lattice  hydrogen  molecules. 


Figure  5 

The  consequences  of  the  structural  details  provided  above  are  manifested  in  the  average  total 
impurity-pI^  interaction  energy.  These  values  are  listed  in  Table  1 .  In  all  the  cases  examined  here, 
the  orientation  dependent  potential  yields  significantly  lower  (more  attractive)  average  potential 
energies.  The  explanation  for  this  remains  that  orientation  dependent  Al-pHj  potential  induces 
some  distortion  of  the  nearest  neighbor  ligand  shell  to  allow  the  A1  p-orbital  to  orient  itself  to 
minimize  the  repulsive  interaction. 


The  difference  in  the  average  potential  energy  values  between  the  vacancy  and  the  no-vacancy 
case  indicate  that  the  presence  of  the  vacancy  helps  to  stabilize  the  system,  especially  when  the 
orientation  dependent  potential  is  used.  When  the  spherical  potential  is  used,  the  difference 


System 

<Vai-h2>  (cm-) 
angle  dependent 

<Vai-h2>  (cm-) 
spherically  averaged 

-263.41 

-177.37 

Impurity  w/o  vacancy 

(3.22) 

(1.26) 

-274.86 

-179.52 

Impurity  w/  vacancy 

(2.37) 

(1.65) 

between  the  average  energies  is  less  than  the  statistical  errors,  so  that  the  extent  of  the  stabilization 
due  to  the  vacancy  is  unclear.  In  our  previous  study  with  a  boron  impurity,  it  was  found  that  the 
average  potential  ener^es  for  the  vacancy  and  no-vacancy  case  differed  less  than  the  statistical  error, 
both  when  the  orientation-dependent  and  spherically-averaged  potentials  were  used.  Hence,  it  is 
clear  that  the  when  the  orientation-dependent  potential  is  used  to  describe  the  AI-/7H2  interaction,  the 
embedded  A1  atom  is  further  stabilized  by  the  presence  of  a  nearest  neighbor  vacancy. 
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We  continued  to  work  on  the  determination  of  the  recombination  rate  of  two  A1  atoms  in  the 
solid  hydrogen  matrix.  This  was  very  important,  since  it  allowed  us  to  better  assess  how  good  of  a 
candidate  AI  was  for  a  possible  dopent  for  the  IBEDM,  Constrained  path-integral  molecular 
dyn^cs  simulation  methods  were  used  to  calculated  the  centroid  potential  of  mean  force  by 
varying  the  centroid  of  the  reaction  coordinate  from  a  large  value  up  to  where  the  recombination 
event  occurs.  Path-inte^al  quantum  transition  state  theory  was  then  utilized  to  evaluate  the  rate 
constant  for  the  recombination  of  two  Al  atoms  in  the  solid  hydrogen  matrix. 
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